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ABSTRACT 

Evidence of a new particle with mass ~125 GeV decaying into a pair of tau leptons at 
the Large Hadron Colhder spurs interest in ascertaining its spin in this channel. Here we 
present a comparative study between spin-0 and spin-2 nature of this new particle, using 
spin correlations and decay product directions. The TauSpinner algorithm is used to re- 
weight distributions from gg — )■ 7/Z — )■ r^r~ sample to simulate a spin-2 state exchange. 
The method is based on supplementing the Standard Model matrix elements with those 
arising from presence of a new interaction. Studies with simulated samples demonstrate 
the discrimination power between these spin hypotheses based on data collected at the 
Large Hadron Collider. 
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1 Introduction 



Following the discovery of a Higgs-like state H with mass ~ 125 GeV at the Large Hadron 
Collider (LHC) [H E] , its spin-parity assignment must be examined to verify the true nature 
of this new state. The spin of this newly observed state has recently been discussed [3] 
in the context of its couplings to a pair of vector bosons. However, from an experimental 
point of view, the spin property should be investigated channel by channel, and other 
alternative hypotheses should be investigated and excluded. At the HCP'12 conference, 
the ATLAS fi] and the CMS [5] collaborations reported observed significances of 1.1 a 
and 1.5 0" respectively, for the H — t- r+r^ decay channel. Their corresponding expected 
significances are 1.7 a and 2.5 a, which when added in quadrature are already at the 3 cr 
level. In the present paper, we concentrate on this channel as a possible completion of the 
spin studies. 

Searches for H — > r+r^ decay are challenging because the r neutrino's escape detection. 
Experimental signatures are categorized over multiple channels in terms of observable 
final state decay products. Data from the multi-channel inputs must be compared with 
simulation of large samples of Monte Carlo (MC) events, which includes detector resolution 
and acceptance effects, as well as contributions from background events in the selected 
sample. 

The study of r polarization can provide additional leverage for this search. The 
TauSpinner algorithm [6] provides a mechanism to evaluate the polarization effects of 
T spin. The algorithm based on re-weighting technique can be applied to existing sample 
of simulated MC events, thereby reducing the need for computationally intensive simu- 
lation of independent samples, and has successfully been applied for measurements of r 
polarization in — )■ r^i/ [7] and Z — )• r+r^ [8] decays. 

In present paper we extend the method of TauSpinner by adding contributions from 
new resonances to the amplitude of gg — ?■ 7/Z — )■ r+r^ processes. Our numerical study 
based on exchange of a spin-2 state X is motivated by recent interest in measurement of spin 
properties for the newly discovered Higgs-like particle candidate. In general, contribution 
from other new interactions, such as those arising from an additional Z' boson, can also 
be evaluated in this way. Though the present implementation illustrates re-weighting of 
samples generated with Pythia [9], the method is equally applicable to other MC event 
generators. 

Our paper is organized as follows. In the next section we discuss contributions from 
spin-2 state to the matrix element. Section 3 is devoted to the TauSpinner algorithm and 
the inclusion of new matrix elements into the program. Section 4 is devoted to technical 
tests of the TauSpinner algorithm and stability tests of internal cross-checks. In section 
5, we investigate experimentally discriminating variables sensitive to spin. In section 6, 
we perform a numerical analysis to access the sensitivity to measure the spin properties of 
Higgs-like states. Section 7 presents the Summary and Appendix A closes the paper with 
detailed description of updates to the TauSpinner algorithm. 
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2 Quark level cross section for ^/Z/X production of 
tau pairs. 

In many theoretical models massive objects of spin-2 arise, including KK gravitons [TUlfTT] . 
analogues of the /2 state of QCD in a new strongly-interacting sector p2] or states in four- 
dimensional ghost-free models of massive gravity [13]. For our purposes we will treat the 
spin-2 particle with mass of 125 GeV as a low-energy signature of some unspecified high- 
energy completion of the model. Therefore, we will use an effective Lagrangian formalism 
for a spin-2 field interacting with fermions to calculate the angular distribution in the 
process qq — )• 7/Z/X — >■ r+r" at the lowest level. These quark level calculations are 
implemented in the TauSpinner [6] algorithm, as described in Section 3. 

For a symmetric spin-2 field X^^ with mass Mx, decay width and momentum k the 
propagator reads [Hj as: 

The projector P is given by: 

pM^^,"/3 = + 7^/^/3 ry'^") - Iry'^^T^"/^ + . . . , (2) 

2 3 

where r]^^ is the Minkowski tensor. The terms proportional to the momentum fc, repre- 
sented by dots in the above formula, will vanish when contracted with the on-shell massless 
fermion currents. 

The interactions of with fermions consists of various operators of increasing di- 
mensions, suppressed by powers of some high scale denoted by F. At zero derivative level 
the coupling of X to a fermion current has a form: 

C 3 XI:4j{XlPl + XrPr)^ + h.c, (3) 

where Pl,r = (1 =F75)/2. This form of coupling is similar to an ordinary Yukawa coupling 
of a Standard Model (SM) singlet scalar. Therefore its experimental signatures: angular 
distributions and spin correlations, will be similar to a scalar exchange, with the only 
difference coming from the spin-2 propagator. If the X^'^ were of gravity or extra dimension 
origin, the couplings Aj for light fermions would naturally be suppressed by the fermion 
mass, A ~ m/F. Therefore, in the discussion of the r-lepton pair production via Drell-Yan 
process in pp collisions such couplings will be ignored. 

At dimension 5 level the coupling of X'^'^ to a fermion bilinear is given by [T5| [16] : 

Cb'- l^X^-la'^ijLi^^d, + ^,d^)iJL + l3\dMl^. + ^^.^lv)i^L] + {L ^ R) + h.c. (4) 

4 r 

The other possible dimension 5 coupling ~ X'^ipc^ip will be ignored since for the on-shell 
fermions can be reduced to the form in eq.Q. 
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The couplings a^'^, P^'^ are model dependent. Although we do not attempt to con- 
struct any specific model, we assume that the couplings are quark- and lepton-flavor diag- 
onal and, following Ref. [16], we make a simplifying assumption: 



a 



L,R 



L,R 



(5) 



At tree level the process qq — )■ r+r proceeds via s-channel '-f/Z/X exchanges. The 
angular distribution in the CM frame can be written as a sum of the SM contribution: 



da 



SM 



dcos6' 



3847rs 



I 2 I 2 n 2 n 2 



^ 2 R 2 
+ ( 9^ap 9^rr + 9- 



'•yqq J-/TT 

(M, 



H 2 I 2 

■yqq 9'^/tt 

s) 



1 + COS 6 
(1 -cos^)^ 



1927r (M| - s)2 + M|r| 

{9^qq 9Zqq9^TT9ZTT + 9^qq9Zqq9yTT9z 

I f L L R R I R R L L 
"T \9'yqq 9 Zqq9-yTT9 Ztt 9^qq9 Zqq9^TT9 Ztt 



1 + COS ( 
COS 9 



3847r (M| - sf + M|r| 



I 2 I 2 
9Zqq 9z TT 



R R ^ 
9Zqq 9z TT 



£^2^2 H L ^ 
9Zqq 9 Ztt + 9Zqq 9z TT 



1 -|- COS ( 

cos^^ 



(6) 



and a new term from the X particle exchange, including its interference with the SM 
contribution. For a real field X and real couplings 9xff reads: 
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In the above expressions 6 is the CM scattering angle of r with respect to the incoming 
quark momentum, and S'a// (with i = L,R) denote the chiral couphngs of A = 'y, Z,X 
bosons to a fermion current //. In particular, for the SM couplings we have gl^jjr = Qf, 
g^ff = cosdw ^'^f sin^ 9w), with g weak coupling constant, and Qf, T| are the fermion 
electric charge and third component of weak isospin. 

The angular distribution due to the X state exchange depends on its couplings. For 
general chiral couplings, the diagonal term is a fourth order polynomial in cos 6, and the 
interference a third order polynomial. For vector or axial-vector type couplings (s'x// = 
±g^jjr), the expressions in eq.([7]) simplify and the diagonal term exhibits a well known 
angular distribution 1 — 3 cos^ 6 + 4 cos^ 6, while the interference with the SM contribution 
behaves as cos^ 6 and vanishes after angular integration. This can be observed in Fig. [6] 
discussed later in the paper. For general couplings the interference survives the angular 
integration. 

The user may take his/her preferred scenario and modify the parameters, including the 
SM ones, in these new currents. For the sake of numerical comparisons, we use F = 1000 
GeV to study the effects of the spin-2 state X at 125 GeV. The width of this state is taken 
to be 1.5 GeV. 

Unless otherwise mentioned, for the default comparison, we set both the left-handed 
and right-handed couplings to have the strengths of unity. We also study two alternative 
scenarios of pure left-handed or pure right-handed couplings labeled in subsequent discus- 
sions as L and R respectively. For these alternative models, the non-vanishing coupling 
constants are increased by \/2 to allow direct comparison with the model {L + R) where 
both couplings contribute equally. 

3 Algorithm case of non-Standard Model weight 

Basic functionality of TauSpinner relies on calculation of effective Born cross sections. It 
is documented in detail in Ref. [B]. For the extension of TauSpinner algorithm, starting 
with the case of production of gg — 7/Z — > r+r~ events, we replace the effective Born 
level contribution with the one where addition amplitudes due to the extra interaction 
are added. At present, the algorithm assumes that the resulting new cross-section has 
contributions from polynomials at most of the 4*'* order in cosine of the scattering angle 6. 
Spin correlations between the produced r"*" and r~ remain as in the case of 7/Z production, 
but differ only in angular dependence of the r polarization. Thus, the non-Standard Model 
weight can easily be calculated using the algorithm described in Appendix B.2 of [6]. 
With minor modifications as interfaced with the Tauola++ MC event generator [17], the 
algorithm is described below: 

• Initialize Tauola++ 

• Initialize TauSpinner 

At this step, user provides two additional flags, nonSM2 and nonSMN. The first one 
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activates calculation of non-Standard Model weight. Then effects due to spin cor- 
relations, described by the ratio of WTl and WT^, and effect on angular dependence 
in T^T~ production is introduced with WT3. The factor due to the ratio of angle- 
integrated cross-section of SM and nonSM case is removed from WT3 if nonSMN is set 
to 1. All other aspects of TauSpinner remain unchanged. 

• Read the data files 

• Calculate the weights 

For the new non-Standard Model case, this step is extended: 

— As in default case, calculate spin weight WTl for Standard Model spin correla- 
tions 

— Calculate spin weight WT2 for spin correlations calculated in non-Standard Model 
case. Use setNonSMkey (1) to switch to non-Standard Model calculation mode. 

— Calculate non-Standard Model weight WT3 for effects on cross section; use double 
getWtNonSMO for its calculation. 

— return WT = (WT2/WT1) * WT3 

An example of such a calculation is given in our tau-reweight-text . cxx. Weight 
calculation is prepared for the case when the generated sample has already the SM spin 
correlations taken into account. User is free to use his own version of the function providing 
the Born level r pair production and featuring another assumptions for SM extension. The 
necessary methods are described in Appendix |X] 

4 Technical tests 

Before any use of our program is started, some technical tests must be performed to verify 
if the function used in the implementation of non-Standard Model interaction is proper 
from the point of view of our program requirements. The first is to check the user provided 
non-Standard Model Born cross- sect ion. In particular, it has to be checked that the same 
conventions for input parameters as those used in TauSpinner and Tauola++ are chosen. In 
our case, for the model as described in previous sections, adjustment of signs for spin states 
had to be introduced. The C++ function nonSM_adopt is prepared for the adjustment of 
conventions. 

An arrangement to verify the proper matching is prepared (an overall y/s dependent 
factors cancel out). Special printouts from TAU0LA/TauSpinner/src/spin2 . cxx, are ac- 
tivated in the DEBUG mode of TauSpinner. Comparison of r polarization calculated 
from TauSpinner default implementation and the ones from user prepared inputs are then 
printed. The corresponding DEBUG mode printout looks as follows: 

^WTl is the spin weight of Standard Model and WT2 is when non-Standard Model interactions are 
switched on. WTl is the default spin weight that TauSpinner returns. 
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\fs [GeV] \fs [GeV] 

Figure 1: The difference between the r polarization calculated analytically from formulae 
as described in Section [2] but with new effects switched off, and the default implementation 
in the Tauola++ are plotted for fixed scattering angle cos^ = 0.3 as a function of y/s. The 
plot on the left is for up quarks, and the plot on the right is for down quark. 
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Figure 2: Angular dependence of the r polarization in the rest frame of hard process is 
shown for virtuality fixed at 125 GeV, for up quarks on the left plot and down-quarks on 
the right plot. The black solid line is the SM contribution from default implementation in 
Tauola+-(-, and the red dashed line is the SM contribution but using nonSM package. The 
effect of nonSM interaction is also shown on the same plots by the blue dotted lines. 
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test of nonSM Born nonsm2=0 
ide,svar,costhe=l 292.05 0.74348 
sni=0. 529759 sm (new) =0 . 530185 nsm=0 . 530185 
sm and smn should be essentially equal. 

The sm and smn denote r polarization as obtained from the Standard Model respectively 
by the method of TauSpinner and the one of the user (the third quantity, nsm, provided by 
the user includes additional interaction, which for svar = s = 292.05 GeV^ is consistent 
with zero). 

Graphic representation of the above tests are provided in Fig. [1] as a function of ^/s 
for a fixed scattering angle 6, and in Fig. [2] for fixed -y/s and as a function of 6. In 
Fig. [U consistency up to per mille level on the r polarization is obtained by replacing the 
TauSpinner standard method of Born calculation with the one as defined in Section 2. For 
nonSM2=0 agreement between the user method and TauSpinner default should be obtained 
as in our case. This is necessary to ensure the non Standard model effect can be correctly 
implemented with nonSM2=l. In Fig. [21 distributions from all SM contribution as well 
as new contributions to the r polarization are superimposed on the same figure. Good 
enough agreement demonstrates the validity of this technical test. Further fine tuning of 
parameters and schemes (e.g. constant or running Z width) in Tauola++ is not necessar}!^. 

Table [1] represents further tests of r polarization, where we check that proper spin 
states of r are provided by the TauSpinner algorithm using the getTauSpinO method 
(see Ref. [6] for definition). Average r polarizations are shown when the virtuality of r 
pair is required to lie within a ± 3 GeV window centered around 125 GeV in the top 
row. For the bottom row, a weight = WT - 1 is used instead of the cut on virtuality. 
Nonetheless, the r polarization still includes contributions from SM via the interference 
effects between the SM and non-SM interactions. The following cases are monitored: SM 
contribution (second column), non-SM contributions for the case of (L + R) couplings 
as in Section 2 (third column), only left coupling for non SM (fourth column) and only 
right couplings (fifth column) . The use of the weight = WT - 1 is of interest for study of 
dedicated sub-samples. 



Selection 


Z T+T- 


Z/X{L + R)^ T+T- 


Z/X{L) T+T- 


Z/X{R) T+T- 


125 ± 3 GeV 


-0.448 ± 0.001 


-0.354 ± 0.001 


-0.521 ± 0.001 


-0.071 ± 0.001 


WT - 1 




-0.255 ± 0.001 


-0.574 ± 0.001 


0.130 ± 0.001 



Table 1: Average r polarization for the SM and non-SM contributions as calculated from 
helicity states attributed by standard TauSpinner method getTauSpinO . 

Second test is performed to check the implementation of non-SM effects when depen- 
dence of the cross-section on r angular distribution and spin effects are included, but 

^One has to bear in mind that TauSpinner constructs the hard process kinematic configuration differ- 
ently than Tauola++. See Ref. [6] for discussion of potential but minor systematic errors. 
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Figure 3: Distributions of invariant masses of the r pair are shown for the SM contribution 
in sohd black hne on the left plot. The almost overlapping dotted blue line includes 
the effect from non-SM contribution, but angular integration dependent contributions to 
the cross section are removed from the weight. The right hand side plot visualizes the 
difference, obtained by applying a weight = WT - 1. 




Figure 4: Distributions of invariant masses of the r pair are shown for the SM contribution 
in solid black line on the left plot. The dashed red line includes the effect from non-SM 
contribution. The right hand side plot visualizes the difference, obtained by applying a 
weight = WT - 1. 
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those arising from integrated Born level cross section are removed. In the TauSpinner 
algorithm, we integrate the distribution over cos 6' for the SM cross-sections at the Born 
level and calculate the effect due to non-SM interactions. If the calculations are correct, 
these two contributions from the new matrix element should compensate each other in 
the distribution of r pair invariant mass rrir+T- , as visible in Fig. |3l The black solid his- 
togram describing the SM contribution is taken from Tauola++ by setting both the flags 
nonSM2 and nonSMN to be equal to 0, whereas for the blue dotted line the new current 
from TauSpinner is used but its effect on cross section are removed. It is tested by setting 
both flags, nonSM2 and nonSMN, to be equal to 1. Change of angular distributions resulting 
from the weight, but integrated over cos 6', results in increased statistical fluctuations as 
visible on the right hand plot in Fig. [3] obtained using weight = WT - 1. The amplitude 
of this procedure is consistent with zero up to per mille level, as expected from the size of 
agreement observed in Fig. [H 

The differential effect of not integrating over the angular dependence due to non-SM 
interactions recovers the expected peaking structure around 125 GeV as shown in the red 
dashed histograms in Fig. |U which are obtained by setting the flags nonSM2 and nonSMN 
to be equal to 1 and 0, respectively. 



5 Spin sensitive observables 

The main feature of TauSpinner is that it can work on previously generated and stored data 
files of the simulated events. For our purposes, we study — )• Z — )• r+r~ generated at y/s 
= 8 TeV using PythiaS [H] MC generator, with r decays simulated by Tauola++ [T7]. To ap- 
ply non-Standard Model weight, TauSpinner was usec^l. The distributions for X — )■ r+r^ 
are obtained by re-weighting the corresponding distributions from Z — )■ r'^r^ samples with 
weight = WT - 1. The samples for gluon-gluon fusion and vector-boson fusion production 
of the Higgs are generated using POWHEG-BOX [181 EH] MC generator interfaced to PythiaS 
for showering and hadronization, while PythiaS is used for the vector-boson associated 
Higgs production. The sample of H ^ t^t^ events are obtained by summing these three 
samples weighted by their respective cross- sect ions [201 EI] • 

We study the variable cos(6'*) |22], describing the average scattering angle between 
the observable final state products from r+ and r~ decays, respectively. In the leading 
approximation, it corresponds to scattering angle in the rest-frame of hard process. In the 
laboratory frame, this variable is defined as: 

sin 6*^ cos 6''^''^ -|- sin 6*"^^ cos 6''^ , , 

cosf6' ) = — 



sin 6'^'' sin 6'^ 

where ^"^^ is the angle between the decay products from r"*" and negative z-axis, and 0'^ 
is the angle between the decay products from r~ and positive z-axis. 

Sensitivity to forward-backward asymmetry is enhanced by flipping the sign of the 
variable cos(6'*), whenever the vector sum of momentum of the visible r daughters is 



•^From Tauola++ version vl . 1 . 0, TauSpinner is provided as a part of the Tauola++ distribution. 
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Figure 5: The distributions of cos(6'*) and cos(^*') variables are shown for Z — )■ r~r"^ 
events before any cuts (on the left) and cuts on truth level invariant mass of the r+r" 
system restricted to lie within a ± 3 GeV window centered around 125 GeV (on the right) 
All distributions are normalized to unit area. 
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Figure 6: The distributions of cos(6'*') are shown for Z — )■ r^r^, X — )■ r+r^, and H — )■ 
r^r~ events after cuts on truth level invariant mass of the r+r^ system restricted to lie 
within a ± 3 GeV window centered around 125 GeV (on the left). For the same selection 
criteria, the distributions of cos(^*') are shown (on the right) for X — t- r^r"*" events for three 
choices of coupling constants as described in the text. All distributions are normalized to 
unit area. 
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aligned towards the negative z-axis. This defines the variable cos{6*'), which is sensitive 
to the spin of Z/X/H. 

The distributions of cos(6'*) and cos(6'*') are shown in Fig. [5]for Z — )■ r^r^ events. The 
left plot in Fig. [5] compares the distributions before any cuts. The right plot compares the 
distributions after a cut on virtuality has been applied by selecting events with the truth 
level invariant mass of the r+r~ system restricted to lie within a ± 3 GeV window centered 
around 125 GeV. 

The distributions of cos(^*') for X — )■ t^t^, and H — r"*"r~ events are shown on the 
left plot in Fig. El The effect of spin is clearly visible in this variable which shows striking 
difference as compared to Z — )■ r+r^ events, also shown in the same plot. The virtuality 
of all these processes are chosen to lie within a ± 3 GeV window centered around 125 GeV 
by applying a cut on the truth level invariant mass of the r+r^ system. 

For the same selection criteria, the distribution of X — )■ r^r+ events are shown on the 
right plot in Fig. [HI for different choices of coupling constants. As expected from Table (H 
the asymmetry in this variable is strongest for the case of pure left-handed coupling, and 
weakest for the case of pure right-handed coupling. 

6 Test for spin of Z/X/H 

A key feature of r polarization in Z — t- t^t~ events is its dependence on virtuality of the 
hard interaction. To select the Z — > r+r~ events as relevant backgrounds for a Higgs-like 
particle with a mass of 125 GeV, we select events inside an appropriate window on the 
collinear mass (m^°"^_) reconstructed from the r decay products [23]. This is sufficient 
for our illustrative purposes, and is comparable to computationally intensive techniques 
developed for di-tau mass reconstruction [21] . 

Our selection criteria are similar to the ones as in Ref . [25l [26] , and those used by the 
ATLAS ^ and CMS [5j collaborations. The following requirements are applied to select 
sample of events enriched with a new resonance produced at 125 GeV: 

• pt for each of the visible daughters of the tau's are required to be greater than 20 
GeV and lie within an acceptance of \ri\ < 2.5, 

• the missing transverse momentum (E^^^^), defined as the of the vector sum of the 
neutrino's momentum, is required to be greater than 20 GeV, 

• the transverse mass of the system comprising of E^^^^ and the visible daughter with 
smaller defined as ■\/2pt-£'t"^^(1 ~ cos A0) is required to be less than 50 GeV, 
where A0 is the angle between the directions of i?™** and visible daughter with 
smaller in the plane perpendicular to the beam direction, 

• the cosine of the 3-dimensional opening angle between the two daughters is required 
to be greater than -0.9, 
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Figure 7: The distributions of cos(^^*') variable are shown for Z — t- r+r~ events with cuts 
on reconstructed level quantity m'^+^_. The plot on the left corresponds to ± 10, ± 5 and 
± 20 GeV windows centered around 125 GeV. The plot on the right corresponds to ± 
10, ± 5 and ± 5 GeV windows centered around 125, 110 and 140 GeV, respectively. All 
distributions are normalized to unit area. 
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Figure 8: The distributions of cos(^*') are shown (on the left) for Z — > r+r", X — i- r+r~, 
and H — )■ t~^t~ events after cuts as described in the text. The distributions of cos(^*') are 
shown (on the right) after cuts for X — )■ t~^t~ events with different options for coupling 
constants as described in the text. All distributions are normalized to unit area. 
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• the difference in azimuthal angles between the two daughters is required to be less 
than 2.9, and 

• m'^+^_ is required to lie within a ± 10 GeV window centered around 125 GeV. 

The stability of the distribution of the discriminating variable cos(^*') for Z — )■ r+r^ 
events with respect to selected region of virtuahty is studied in Fig. [71 The requirements 
on the window on m'^J^^_ is varied, keeping all the other above-mentioned selection require- 
ments the same. The left plot corresponds to ± 10, ± 5 and ± 20 GeV windows centered 
around 125 GeV. The right plot corresponds to ± 10, ± 5 and ± 10 GeV windows cen- 
tered around 125, 110 and 140 GeV, respectively. For the case when virtuality is centered 
around 125 GeV, the distributions are compatible, while the shape is altered with shift of 
the center of the window. Thus, these differences in shapes can be taken as estimates of 
systematic uncertainties on the modeling of the background from Z — t- r^r+ events. 

The left plot in Fig. [8] shows the distribution of cos(^*') after applying these selection 
requirements to sample of Z — )■ t~^t~ , X — )■ t~^t~, and H — )■ r^r~ events. For illustrative 
purposes, the distributions of cos(6'*') are shown in the right plot of Fig. |H]for X — t~^t~ 
events passing these selection criteria corresponding to the three different choices of the 
coupling constants. 

Visible separation between the shapes of if — )■ r+r~ and X — )■ r+r" events are ob- 
served, which is similar in size for all the three choices of coupling constants studied. The 
distributions from H — > r+r~ events are incompatible at the 99% confidence level with 
those arising from X — )■ r+r^ events for a sample of 500 Higgs-like events, which corre- 
sponds to roughly twice the number of signal events observed with 13 fb^^ at pp collision 
energy of 8 TeV by the ATLAS collaboration However, the separation power depends 
on the presence of backgrounds and choice of selection criteria, which alters relative effi- 
ciencies of signal and background events. Further discrimination power can be obtained 
by categorizing the r+r^ decay modes depending upon the observed final state particles. 

7 Summary 

In this note we have presented a new development of the TauSpinner that is capable 
of addressing possible extension of the SM based on new physics model that results in 
contribution from a spin 2 state X. To simplify estimation of observability of such a state 
at the LHC, extension for the TauSpinner algorithm to manipulate r pair final states in 
previously generated MC samples is proposed which employs the method of event weights. 
Weights calculated with TauSpinner feature an implementation of amplitudes from beyond 
the Standard Model. Not only spin correlation effects as in the previous versions, but also 
effects modifying the angular distributions of lepton directions and size of the cross- 
section can be studied in this way. 

An example of the installation of our new matrix elements into the TauSpinner algo- 
rithm is presented. Distributions validating correctness of the installation are discussed 
and are found correct. 
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We then study distributions of experimentally observable quantities sensitive to the spin 
of Z/X/H. In semi-realistic conditions, we study the impact of the new interaction and 
its signature. We found the approximate hard process scattering angle 9*' reconstructed 
from observable directions of r decay products useful for that purpose. Following the 
experimental searches performed with half the dataset expected to be collected at \/s = 
8 TeV in 2012, we apply selection criteria that enrich the sample of events with a Higgs- 
like state with mass ~ 125 GeV. With double the number of Higgs-like events as can be 
expected using full dataset in 2012, we find distributions that are incompatible between 
H — )• r+r" and X — )■ t~^t~ events at the 99% confidence level. Detailed experimental 
study with all background contributions, which include detector resolution and acceptance 
effects, as categorized over the different r^r^ decay modes is expected to improve the 
discrimination power. 

Robustness of the method was demonstrated here, and the first results are of potential 
interest. Method is also straightforward to extend to other cases such as Z' etc. Alter- 
native production mechanism of spin-2 states, eg. via gluon-gluon fusion, may also be 
implemented by appropriate re-weighting of the angular dependence from corresponding 
matrix elements. A key aspect of our TauSpinner algorithm is that computationally ex- 
pensive simulation of independent MC samples is not necessary for study of new physics 
interactions. 

A TauSpinner - changes introduced in version 1.2 

Since its first public version, described in [6], two new updates to TauSpinner has been 
introduced. First, TauSpinner has been merged into Tauola++[l7j distribution and now, 
while working on this paper, it has been extended to add new functionality. In this section 
we list the changes between version 1.0 and 1.2. 

• Merging with Tauola++ 

TauSpinner now comes as an additional library to Tauola++ distribution. Tauola++ 
configuration scripts have been updated to accomodate this setup. As of writ- 
ing this paper, Tauola++ vl . 1 . 1, featuring TauSpinner vl.2 has been installed 
in GENSER^ |28] database. 

• Two new initialization options - nonSM2 and nonSMN 

The nonSM2 flag turns on non-Standard Model weight calculation. The nonSMN flag, 
combined with nonSM2, allows for calculation of corrections to shape only. 

• New functions added. 

An example examples/tau-reweight-test . cxx has been updated to show how func- 
tions described below can be used in case of spin-2 calculation described in this paper. 

— set_nonSM_born( double (*fun) (int ID, double S, double cost, int HI, 
int H2, int key) ) 
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Sets function for user-defined born, including new physics. The parameters of 
this new function are described in include/TauSpinner/nonSM.h as well as in 
the example program. 

— void setNonSMkeyCint key) 

Sets the value of nonSM2 flag. Allows turning non-Standard Model calculation 
on and off for comparison between different models. 

- double getWtNonSMO 

Returns non-Standard Model weight WT3 calculated for the last event processed 
by TauSpinner. 

In our example examples/tau-reweight-test . cxx fortran function is provided to cal- 
culate quark level Born cross section where new physics effects can be switched on and off. 
Physics model described in the previous sections at the level of quark level annihilation 
into pair of tau leptons is used. Function: 

REAL*8 FUNCTION D1STJWK(1D,S,T,H1 ,H2,KEY) 

is used in our program with the help of the C++ function: 

double nonSM_adopt (int ID, double S, double cost, int HI, int H2, int key) 

Its use is initialized with the method set_nonSM_born( nonSM_adopt ) ; 
Other user defined function can be used in the same way. 
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